function a = aitken( a0, a1, a2 )
% function a = aitken( a0, a1, a2 )
%   arguments:
%     a0, a1, a2	= three successive numbers in the sequence
%   returns:
%     a				= the aitken estimate 
	if ( abs(a2 - 2*a1 + a0) < 1e-9 && ( a2 - a1 )^2 < 1e-10 )
		a = a2 - ( a2 - a1 )^2 / ( a0 - a2 );
	elseif ( abs(a2 - 2*a1 + a0) < 1e-9 )
		a = a2;
	else
		%a = a2 - ( a2 - a1 )^2 ./ ( a2 - 2*a1 + a0 );
        a = ( a2*a0 - a1^2 ) / (a2 -2*a1+a0);
 	end
end

